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ABSTRACT 

Context. Early B-type stars are invaluable indicators for elemental abundances of their birth environments. In contrast to the sur- 
rounding neutral interstellar matter (ISM) and H ii regions their chemical composition is unaffected by depletion onto dust grains and 
by the derivation of different abundances from recombination and collisional lines. In combination with ISM or nebular gas-phase 
abundances they facilitate the otherwise inaccessible dust-phase composition to be constrained. 

Aims. Precise abundances of C, N, Mg, Ne, Fe in early B-type stars in the Orion star-forming region are determined in order to: a) 
review previous determinations using a self-consistent quantitative spectral analysis based on modem stellar atmospheres and recently 
updated model atoms, b) complement results found in Paper I for oxygen and silicon, c) establish an accurate and reliable set of stellar 
metal abundances to constrain the dust-phase composition of the Orion H ii region in Paper II of the series. 

Methods. A detailed, self-consistent spectroscopic study of a sample of 13 narrow-lined B0V-B2 V stars in OriOBl is performed. 
High-quality spectra obtained with FIES@NOT are analysed using both a hybrid non-LTE method (classical line-blanketed LTE 
model atmospheres and non-LTE line formation) and line-profile fitting techniques, validating the approach by comparison with re- 
sults obtained in Paper I using line-blanketed non-LTE model atmospheres and a curve-of-growth analysis. 

Results. The comparison of the two independent analysis strategies gives consistent results for basic stellar parameters and abun- 
dances of oxygen and silicon. The extended analysis to C, N, Mg, Ne, and Fe finds a high degree of chemical homogeneity, with the 
Icr-scatter typically adopting values of 0.03-0.07 dex around the mean for the various elements. Present-day abundances from B-type 
stars in OriOBl are compatible at similar precision with cosmic abundance standard values as recently established from early-type 
stars in the solar neighbourhood and also with the Sun. 



Key words. Stars: abundances 
Associations 



Stars: atmospheres — Stars: early-type — Stars: fundamental parameters — Open clusters 



> 

in 
m 

o 



X 



1. Introduction 

The Orion complex, containing the Orion molecular cloud and 
the Orion OBI (OriOBl) association, is one of the nearest 
Galactic regions with ongoing star formation, located at a dis- 
tance of about 400 pc. It is unique for studying young massive 
stars, interstellar matter, and their interaction via feedback pro- 
cesses. The Orion complex can be regarded as the best reference 
for any further study of other star-forming regions. Its proxim- 
ity and brightness have facilitated many detailed investigations, 
that yielded important discoveries and a wealth of information 
on many astrophysical phenomena over the past decades. 

iBlaauwl d 19641) divided Ori OB 1 into four subgroups of stars, 
namely la, lb, Ic and Id, distinguish ed by their differen t lo- 
cation in the sky and by their ages. iBrown et al.l (11994 de- 
rived mean ages of 11.4+1.9, 1.7+1.11, 4.6+2, and <lMyr 
for subgroups la to Id, respectively. The youngest subgroup 
OriOBl Id is associated with the Orion nebula (M42), the 
most studied H ii region and the closest ionized nebula to the 



* Based on observations made with the Nordic Optical Telescope, op- 
erated on the island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astroflsica de Canarias. 
** Visiting researcher at Dr Karl Remeis-Stemwarte & ECAP, 
Universitat Erlangen-Niimberg. 



Sun. The correlation between the ages of the stellar subgroups, 
their location, and the large scale structures in the interstellar 
medium around Orion OBI have been interpreted as features 
of seq uential star formation and the impact of type-II super- 
novae dRevnolds & Ogdenll979l:ICowie et alJ1979tlBrown et al.l 
11994 . 

Early B-type stars in OriOBl are important astrophysical 
objects that can be used to obtain present-day chemical abun- 
dances in terms of galactic chemical evolution. Moreover, the 
possibility of comparing the derived stellar abundances with 
those resulting from the spectroscopic analysis of the Orion neb- 
ula makes the tandem massive stars in Ori OBI and M42 a key- 
stone in the investigation of the reliability of B-type stars and H ii 
regions as tracers of present-day abundances in the Universe. 

First abundance deteriidnations in B-type stars in Ori OBI in 
the 90;s (iGies & Lamber"ill992HCunha & Lamberjl 199211 1994 



iKilianll 19921: iGummersbach et al.lll998h indicated a large scatter 
up to 0.4-0.5 dex for O and Si, but also > 0.2 dex for other el- 
ements such as C, N, an d Fe. At that time, the most ex tensive 
work was carried out by ICunha & Lambert] (Il992[ Il994l) . here- 
after CL92 and CL94, respectively. They derived C, N, O, Si, 
and Fe abundances of 18 B-type main sequence stars from the 
four subgroups comprising the OriOBl association. They ob- 
tained a similar scatter in O and Si abundances as in the other 
cited works (based on smaller samples of stars), but also found 
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that the highest abundances were associated to some of the stars 
in the youngest (Id and some Ic) subgroup^ This abundance 
pattern was interpreted as an observational evidence supporting 
a scenario of induced star formation in which the new genera- 
tion of stars are formed from interstellar material contaminated 
by type-II supemovae (SNII) ejecta. 

This picture has changed substantially after new studies of 
ox ygen and silicon abundanc es of early B-type stars in Ori OB 1 
by ISimon-Diaz et all ( l2006h and ISimon-Diaj ( 1201 Ol hereafter 
Paper I). Both stellar parameters and chemical abundances were 
shown to have different values than found previously, for the 
same stars. As a consequence, high chemical homogeneity of O 
and Si as derived from early B-type stars the Orion star-forming 
region is implied, and no clear indication for contamination with 
SN II nucleosynthesis products is found. 

The results from Paper I for the mean abundance and stan- 
dard deviation for O and Si are in full agreement with recent 
findings of homogeneous present-day chemical abundances in 
the solar neighbourhood (Przvbilla, Nieva & Butler 2008, here- 
after PNB08). There, a small but representative sample of early 
B-type stars was analysed, obtaining a significantly smaller scat- 
ter for C, N, O, Si, Mg, and Fe abundances than commonly found 
by previous studies. 

Some reasons causing improvements in c hemical abundance 
determi nations of B-type stars can be found in lNieva & Przvbillal 
(l20T0allbl) . and references therein. In brief, these are a conse- 
quence of the minimization of systematic uncertainties of obser- 
vational origin and from the spectral modeling and spectral anal- 
ysis, with two main contributions: (a) the application of a self- 
consistent quantitative spectral analysis, i.e. to derive stellar pa- 
rameters and chemical abundances from the spectrum only, with- 
out resorting to any extra ingredient like a photometric effective 
temperature estimation, and (b) the reliability of the atomic mod- 
els and the selection of appropriate diagnostic lines for the stellar 
parameter and abundance determination. 

In view of the results obtained in Paper I and PNB08 a re- 
vision of abundances of other important elements in Ori OBI 
stars is warranted. The main aim of this third paper of the se- 
ries, in combination with Paper I, is to establish an accurate and 
reliable set of abundances for C, N, O, Ne, Mg, Si, and Fe in 
early B-type stars in the Orion star-forming region, providing 
a basis to perform a thorough comparison between stellar and 
nebular chemical abundances i n the Orion star-forming region 
ISimon-Dfaz & Stasihskal (|20I II hereafter Paper II of this series). 

The analysis presented in this paper is based on recently up- 
dated model atoms for non-LTE calculations (see Table[T]i which 
have been already applied in a small sample of early B-type 
stars in the solar neighbourhood (PNB08) and more recently 
in a larger sample (Nieva & Przybilla, in prep.). These model 
atoms are not available at present in the stellar atmosphere code 
used in Paper I, therefore, we decided to base the analysis on the 
codes and techniques described in PNB08. For consistency, we 
performed the whole analysis from scratch (also including the 
stellar parameters, and O and Si abundance determination). The 
comparison of results referent to stellar parameters, O and Si 
abundances obtained by two independent analyses will be very 
valuable to reinforce the conclusions presented in Paper I and/or 
investigate possible biases originated by the use of diff'erent stel- 
lar atmosphere codes and methodologies. 



' For other elements the scatter was also relatively large (e.g. 0.2, 
0.45, and 0.35 dex for C, N, and Fe, respectively), but they did not find 
significant variations of abundances across the subgroups. 



The paper is structured as follows: Sect. |2]briefly describes 
the star sample and the observations. The spectrum synthesis in 
non-LTE and the codes used in the present work are described in 
Sect. [3] The self-consistent quantitative spectral analysis is sum- 
marized in Sect. |4] Sect, presumes the results and the compari- 
son with other representative works in the literature. A summary 
is given and the conclusions are discussed in Sect. |6] An exam- 
ple of a global spectrum synthesis for one star of the sample is 
given in AppendixlAl 

2. The star sample 

The star sample consists of 13 early B-type main sequence stars 
of spectral classes BO V to B2 V and low projected rotational 
velocities (vsin!<60 kms ') in the Ori OBI association. The 
spectroscopic dataset used in this study is the same presented 
in Paper I. It refers to high resolution and high signal-to-noise 
ratio spectra obtained with the Fibre-fed Echelle Spectrograph 
(FIES) attached at the 2.5m Nordic Optical Telescope (NOT) in 
El Roque de los Muchachos observatory on La Palma (Canary 
Islands, Spain) in November, 2008. Details of the data reduction 
process and the quality of the spectra can be found in Paper I. 

3. Spectrum synthesis in non-LTE 

The non-LTE line-formation computations follow the methodol- 
ogy discussed in detail in our p revious studies for H, He and C 
(Nieva & Przybilla 20071 l2008h and for N, O, Ne, Mg, Si and 
Fe (PNB08). In brief, a hybrid non-LTE approach is employed, 
i.e. combining classical line-blanketed LTE model atmospheres 
and non-LTE line-formation calculations. This technique pro- 
vides an efficient way to compute realistic synthetic spectra in 
all cases where the atmospheric structure is close to LTE, as it is 
the case for the early B-t ype main sequence stars analysed here 
(iNieva & Przvbillall2007l) . 

The mod el atmospheres were computed with the 
Atlas9 code dKuruc i Il993bh which assumes plane-parallel 
geometry, hydrostatic, radiative and local thermodynamic 
equilibrium, and c h emica l homogeneity. Solar abundances of 
iGrevesse & Sauval (119981) were adopted in the atmospheric 
structure^ computations. Line blanketing was r ealized here by 
means of opacity distribution f unctions ODFs (Kurucz 1993ab 
with solar abundances scaled to lGrevesse & Sauval (.1998.) . The 
model atmospheres were held fixed in the subsequent non-LTE 
calculations. 

Non-LTE level populations and model spect ra were obtained 
with recent versions of Detail and Surface (lGiddingsl[T98l 
iButler & GiddingsllI985h : both updated by K. Butler The cou- 
pled radiative transfer and statistical equilibrium equations were 
solved with Detail, employing an accelerated lambda iteration 
scheme of Rybicki & Hummer (1991). This allows even com- 
plex ions to be treated in a realistic and efficient way. Synthetic 
spectra were then calculated with Surface, using refined line- 
broadening theories. The sources of model atoms adopted for 
the present work are listed in Table [1] Details on the updates are 
given in by Nieva & Przybilla (in prep.). 

Continuous opacities due to hydrogen and helium were con- 
sidered in non-LTE and metal line blocking was accounted for 
in LTE via Kurucz' ODFs. Microturbulence was considered in 
a consistent way throughout all computation steps: selection of 



^ The effect of the 'new' solar abundances (Asplund et al. 2009) on 
the atmospheric stmcture and on the ODF's is weak. 
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Table 1. Model atoms for non-LTE calculations. 



Ion Model atom 

H Przvbilla & Bu tler (2004) 

He i/ii ^zybilla (2005) 

C ii-iv Nieva & Przvbilla (2006. 2008) 

Nn Przvbilla & Butler (2001)* 

O i/ii Przvbilla et al. (2000), B ecker & Butler (1988a' 

Ne i/n Morel & Butler (2008)* 

Mg II Przvbilla et al. ( 200 1) 

Si m/iv Becker & Butler ( 1 990a) 

Fe ii/iii Becker ( 1998), Morel et al. (2007)' 

* Updated by N. Przybilla (priv. comm.) 



appropriate ODFs for line blanketing and line blocking, atmo- 
spheric structure calculations, determination of non-LTE popu- 
lations and formal solution. 

Grids of synthetic spectra for H/He, C, N, O, Mg and Si have 
been computed in order to speed up the analysis. They cover ef- 
fective temperatures from 15 000 to 35 000 K in steps of 1000 K, 
surface gravities from 3.0 to 4.5 in steps of 0.1 dex, microtur- 
bulences from to 8kms ' in steps of 2kms"' (for H/He, Si 
and O microturbulence reaches a value of 12kms ') and metal 
abundances within 1 dex centered on the Cosmic Abundance 
Standard proposed in PNB08 in steps of 0. 1 dex. Hydrogen and 
helium abundances adopt PNB08 values. The lower limit of the 
surface gravity for each value of temperature/microturbulence 
was constrained by the convergence of Atlas9. All grids have 
been successfully tested by reproducing results from PNB08 
within a broad temperature range. Micro-grids varying abun- 
dances only, were computed per star for neon and iron once 
all stellar parameters have been determined with the larger pre- 
computed grids. 

4. Self-consistent quantitative spectral analysis 

The analysis method used in this paper is based on a spectral 
line-fitting procedure^. We aim at deriving atmospheric param- 
eters and chemical abundances self-consistently by reproduc- 
ing several independent spectroscopic indicators simultaneously. 
The atmospheric parameters primarily derived here are the ef- 
fective temperature Tgif, surface gravity \ogg, microturbulence 
macroturbulence ^, projected rotational velocity v sin /, ele- 
mental abundances e(X). 

The quantitative analysis is mostly automatized, but it also 
allows to work interactively in some decisive steps. The inter- 
polation in stellar parameters via spectral line fitting using 
minimization within the precomputed grids is performed with 
the program SpaA The code allows flexible selection of spec- 
troscopic indicators for parameter determinations which may 
vary from star to star upon availability of specific spectral lines. 
Linelists can also be chosen interactively. This is crucial since 
then problematic spectral lines, i.e. too weak, blended, with low 
S/N, in a region where the continuum normalization is difficult 
or with a particular shortcoming in the spectral synthesis can 
be excluded. The local continuum can also be adjusted interac- 
tively. This might slow down the analysis process in cases where 
the linelist is long but the gain in precision is high. 

' Note that we used here a different approach than in Paper I. More 
details can be found in Sect. 15. II 

Spectrum Plotting and Analysing Suite (SPAS. lHirschll2009h . 



Table 2. Spectroscopic indicators for the Tf-g and \ogg determi- 
nation. The boxes indicate ionization equilibria that have been 
quantitatively matched in the analysis. 



# Object H HeiHeii CiiCni Civ Oi On NeiNeii SiiiiSiiv FeiiFeiii 

1 HD36512 » I » » I I « » » I I I I » » I r~ 

2 HD37020 • « » » » | • | » » | » » • 

3 HD36960 • » » » » • » » « » • 

4 HD37042 • « » » » • » » « » • 

5 HD36591 • » » » » | » » | » » « » • 

6 HD36959 • _* »_ • • • • • • • • • 

7 HD37744 • • » » » » » » « » | » » 

8 HD35299 • • » » » » • » » » » 

9 HD36285 • • j *_ j *_ • _* »_ _* »_ 

10 HD35039 • • » » » » • » » « » 

11 HD36629 • • ~i T ~i T • • ~i T 

12 HD36430 • • » » » » • • » » 

13 HD35912 • • 1 T ~ T . • ~ T 



The strategy used for the analysis of our star sample is as 
follows. A first estimation of Teif and log g is obtained by means 
of a simultaneous fit to most H and He lines (an accuracy bet- 
ter than ~5% and 0.1-0.2 dex, respectively, is expected for the 
case of this high quality set stellar spectra). Then the proce- 
dure starts to consider lines of other elements to end up with 
a self-consistent solution for atmospheric parameters and chem- 
ical abundances, also including an estimation of their internal 
errors. In this second step, every element is analysed indepen- 
dently and some interactive iterations to fine-tune the stellar pa- 
rameter determination are needed in order to find a unique solu- 
tion for all indicators. 

Numerous metal lines were analysed simultaneously per 
star/element. The high S/N and high resolution spectra, the op- 
timum continuum normalization, the ample wavelength range 
coverage, the presence of many unblended sharp spectral lines 
and the adoption of our updated model atoms for spectral syn- 
thesis calculations allow us to analyse more lines than in stan- 
dard works. The linelist varies slightly from star to star because 
the selection of good lines to be analysed depends on the tem- 
perature, the quality of the spectra and the wavelength coverage. 
In particular, the line lists used to derive final values of elemen- 
tal abundances has been chosen for each star, in a way that the 
abundance from any line does not differ by more than ~2-3o- 
from average. 

Ionization equilibria, i.e. the requirement that lines from dif- 
ferent ions of an element have to indicate the same chemical 
abundance, along with the fitting of the wings of the H Balmer 
lines facilitate a fine-tuning of the previously derived parame- 
ters. Whenever possible, we use multiple ionization equilibria 
for the parameter determination as the redundancy of informa- 
tion helps to minimize systematic errors. In particular, we have 
found that the use of more than one ionization equilibrium helps 
to better constraint the surface gravity compared to the case 
when only H lines are used. 

The selection of ionization equilibria to be analysed depends 
on the temperature range of the star. Table |2] summarizes the 
available indicator^ for the star sample. Note that the use of 
more than one ionization equilibrium has been always possi- 
ble. In addition, the availability of the Fe ii/iii ionization equi- 
librium allow to also impose a constraint to the metallicity, 
[M/H]:^[Fe/H]. 

' Note that we analyse Siiii/iv lines only because the model atom 
employed here does not produce reliable Si ii lines (see also Paper I). 
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Table 3. Stellar parameters and elemental abundances of the program stars. 







Sub- 


FW 






ADS 












lo 


gX/H+ 12 (NLTE, ADS) 








Target 


SpT 


Gr. 








log 


S 






vsini 


(T 


C 


N 


Ne 


Mo 


Fe 


# 








(K) 


(cgs) 


(K) 


(cgs) 






(km s"' 


) 


















36512 


BOV 


Ic 


DJ /UU 




33400 + 300 


4.30 + 


0.07 


4 + 


1 


20+1 


10 + 2 


8.36 + 0.15 


7.83 + 


0.15 


o 1 A _i_ n 1 n 
o. 10 + u. lU 


/.jl 


7.55 


+ 0.02 


1 


37020 


B0.5V 


Id 




4.Z 


30700 ± 300 


4.30 + 


0.08 


+ 


1 


45 + 3 


20 + 5 


8.40 + 0.07 


7.85 + 


0.09 


no _i_ n no 
o.Uy + U.Uo 


/.o4 


7.53 


+ 0.08 


2 


36960 


B0.5V 


Ic 


zoyuu 


"2 Q 

j.y 


29000 + 350 


4.10 + 


0.08 


4 + 


1 


28 + 2 


20 + 6 


8.39 + 0.07 


7.78 + 


0.08 


o.Vd + U.UV 


/.ol 


7.53 


+ 0.11 


3 


37042 


B0.7V 


Id 


Zy /UU 




29300 + 300 


4.30 + 


0.09 


2 + 


1 


30 + 2 


10 + 3 


8.33+0.11 


8.04 + 


0.08 


0.1 J + u.uy 


n A/1 
/.o4 


7.54 


+ 0.09 


4 


36591 


Bl V 


lb 


27200 


4.1 


27000 + 250 


4.12 + 


0.05 


2 + 


1 


12+1 




8.30 + 0.10 


7.81 + 


0.09 


8.10 + 0.07 


7.61 +0.02 


7.49 


+ 0.12 


5 


36959 


Bl V 


Ic 


25900 


4.2 


26100 + 200 


4.25 + 


0.07 


+ 


1 


12+1 


5 + 1 


8.37 + 0.11 


7.85 + 


0.09 


8.18 + 0.10 


7.66 + 0.04 


7.53 


+ 0.11 


6 


37744 


B1.5V 


lb 


23800 


4.1 


24000 ±400 


4.10 + 


0.10 


+ 


1 


33 + 2 


15 + 2 


8.32 + 0.07 


7.81 + 


0.10 


8.12 + 0.08 


7.58 


7.52 


+ 0.07 


7 


35299 


B1.5V 


la 


23700 


4.2 


24000 + 200 


4.20 + 


0.08 


+ 


1 


8+1 




8.37 + 0.07 


7.81 + 


0.07 


8.12 + 0.09 


7.61+0.07 


7.49 


+ 0.10 


8 


36285 


B2V 


Ic 


20600 


4.0 


21700 + 300 


4.25 + 


0.08 


+ 


1 


11 + 1 


8 + 1 


8.32 + 0.07 


7.77 + 


0.09 


8.04 + 0.08 


7.52 + 0.08 


7.54 


+ 0.10 


9 


35039 


B2V 


la 


19800 


3.7 


19600 + 200 


3.56 + 


0.07 


4 + 


1 


12+1 


7 + 1 


8.34 + 0.10 


7.77 + 


0.09 


8.04 + 0.08 


7.49 + 0.06 


7.48 


+ 0.10 


10 


36629 


B2V 


Ic 


20000 


4.1 


20300 + 400 


4.15 + 


0.10 


2 + 


1 


10+1 


5 + 1 


8.36 + 0.07 


7.82 + 


0.08 


8.02 + 0.09 


7.50 + 0.02 


7.52 


+ 0.11 


11 


36430 


B2V 


Ic 


18600 


4.1 


19300 + 200 


4.14 + 


0.05 


+ 


1 


20 + 2 


10 + 2 


8.32 + 0.08 


7.76 + 


0.10 


8.03 + 0.12 


7.54 + 0.05 


7.50 


+ 0.07 


12 


35912 


B2V 


la 


18500 


4.0 


19000 + 300 


4.00 + 


0.10 


2 + 


1 


15 + 1 


8 + 1 


8.33+0.09 


7.76 + 


0.07 


8.05 + 0.11 


7.50 + 0.05 


7.52 


+ 0.08 


13 



Notes. Target id are HD numbers. Parameters derived with Fastwind (FW, Paper I) and Atlas+Detail+Surface (ADS, this work) are listed. 
Uncertainties in parameters from Paper I are similar than in the present work. 



Table 4. Average metal abundances. 





This work 


Paper I 


PNBP8 


CL94NLTE 


CL94™ 


Sun" 


c 


8.35+0.03 [0.09] 




8.32+0.03 


8.40+0.11 


8.35+0.06 


8.47+0.05 


N 


7.82+0.07 [0.09] 




7.76+0.05 


7.76+0.13 


7.79+0.11 


7.87+0.05 


O 


8.77+0.03 [0.10] 


8.73+0.04 


8.76+0.03 


8.72+0.13 


8.73+0.13 


8.73+0.05 


Ne 


8.09+0.05 [0.09] 




8.08+0.03 


8.11+0.04* 




7.97+0.10 


Mg 


7.57+0.06 [0.05] 




7.56+0.05 






7.64+0.04 


Si 


7.50+0.06 [0.06] 


7.51+0.03 


7.50+0.02 


7.14+0.13 


7.36+0.14 


7.55+0.03 


Fe 


7.52+0.02 [0.09] 




7.44+0.04 




7.47+0.10 


7.54+0.04 



Notes. Average metal abundances, 1-cr uncertainties and intrinsic uncertainties for each element (in brackets) derived in this work are compared 
to values from Paper I, the present-day Cosmic Abundance Standard by PNB08, abundances in NLTE and LTE of 17 stars in Ori OBI by CL94, 
and solar abundances. Protosolar values by AGSS09 (Table 5 in the original paper). Value from 16 stars, from CHL06. 



The fine-tuning of microturbulence, stellar parameters and 
abundances is performed at the same time. In particular, micro- 
turbulence is determined by imposing a zero correlation between 
individual line abundances for some ions. 

The strength of this analysis method is that most im- 
portant sources of systematic uncertainties can be identi- 
fied and_Jhe_jemaining_OTors reduced to minimum values. 
See iNieva & Przvbillal (l2010 i^l5) for more detailed discussions 
on the methodology and systematic errors in the spectral analy- 
sis of this kind of stars. 



5. Our results for the Orion B-type stars and 
comparison with previous viorV. 

Table |3] summarizes atmospheric parameters and elemental 
abundances of C, N, Ne, Mg and Fe derived in this study for 
the 13 sample stars. There, also spectral type and subgroup to 
which the stars belong within Ori OBI are given. The identifica- 
tion numbering of the stars (last column) follows the one used in 
Paper I. The corresponding mean values and 1-cr uncertainties 
derived from the whole star sample are presented in first column 
of Table |4] Note that mean abundances for O and Si obtained 
in this work are also indicated. These abundances (though not 
presented in Table |3]l were derived for completeness, and to be 
compared with those obtained in Paper I (see Sect. l5.lt . 

Table |4] also contains information about results obtained by 
previous studies o f early B-type stars in Ori OBI from Paper I, 
CL92, CL94 and ICunha et all (l2006h . hereafter CHL06, early 



B-type stars in the solar neighbourhood from PNB08, and th e 
set of protosolar abundances proposed bv lAsplund et al.l ( |2009|) . 
hereafter AGSS09. Values from CL94 in TableS]have been com- 
puted from their total sample and are slightly diff'erent than those 
values indicated by them, calculated from a sub-sample of their 
work. A detailed discussion on the comparison of results from 
the analysis performed here and the various cited studies is pre- 
sented in the following sub-sections. 

^From inspection of Tables[3]and|4]it can be concluded that 
the large scatter in abundances found for C, N, and Fe by pre- 
vious works (see Sect. [U has been significantly reduced. In ad- 
dition, the high degree of homogeneity of derived abundances 
for all the analysed elements is remarkable (the 1-cr dispersion 
associated to the mean value of the whole sample is always sig- 
nificantly smaller than the intrinsic uncertainties from the indi- 
vidual analysis, obtained by averaging the individual errors for 
each star listed in Table |3] These values are indicated in brack- 
ets in Table nil. To better illustrate these conclusion, we plot in 
Fig.|2] the individual abundances for the analysed Ori OBI stars 
(black filled squares) together with the 1-cr dispersion of abun- 
dances derived from the whole sample (grey-shaded area), and 
compare them with other studies listed in Table |4] 

5.1. B stars in OhOB1: Paper I 

The same set of observed spectra was analysed in Paper I 
and in the present work. The method used here based 
on Atlas+Detail-i-Surface and the one adopted in Paper I 
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Fig. 1. Comparison of effective temperature and surface gravities (left panel) and oxygen and silicon abundances (right panel) 
derived in Paper I with Fastwind and the present work with Atlas+Detail+Surface. Evolutionary tracks (corresponding to Z=0.02) 
are extracted from lSchaller et al1(ll992l) . Numbers in the left panel identify the stars, as in Table|3]The blue cross in the right panel 
indicates the typical uncertainties in O and Si abundances from individual analyses. 



(Fastwind) employ the same philosophy: the simultaneous de- 
termination of stellar parameters and chemical abundances from 
the spectrum only. However, the codes and the spectral analy- 
sis are different. Results of Paper I are ba sed on the unified at- 
mosphere code Fastwind (IPuIs et al.ll2005h that models the stel- 
lar atmosphere and wind in non-LTE. The input atomic data for 
the line formation of hydrogen, helium, silicon and oxygen are 
not identical to the ones used here. Furthermore, the determina- 
tion of chemical abundances in Paper I was based on equivalent 
widths and the curves of growth, while here direct spectral line 
fits are performed. Finally, the spectroscopic diagnostics used 
for the stellar parameters determination are different: on the one 
hand. Paper I relies exclusively on the H Balmer lines, along 
with the He i/ii and/or Si ii/iii/iv ionization equilibria; on the other 
hand, as described in Sect|4l the analysis presented in this paper 
uses all available ionization equilibria for He, Si, C, N, O, Ne, 
and Fe. 

Stellar parameter and oxygen and silicon abundance deter- 
minations were carried out independently to the work in Paper I, 
without a priori information of final results. The comparison of 
both studies can be found in Tables [3] (stellar parameters) and 
|4](mean abundances and 1-cr dispersion), and Figs. [1] and |2] In 
spite of all differences in the two analyses, the similarity of re- 
sults derived in both works is remarkable. The global agreement 
in the four quantities, Tgft, log g, e(0), and e(Si), is very good 
within the uncertainties. 

Fig.|2]shows a detailed comparison of abundances for the 13 
sample stars with results for O and Si from Paper I (red open 
squares) for the same objects. The analysed stars are separated 
in the various Ori OBI stellar subgroups. No trend in O or Si 
is found, confirming the conclusions of Paper I of chemical ho- 
mogeneity of these elements in the region. We refer the reader 
to Paper I for a similar comparison with results from CL92 and 
CL94. The differences of mean abundances for the sample is 
0.04 dex for oxygen and 0.01 dex for silicon, that within the er- 
rors are negligible, and the 1-cr dispersion of abundances in the 
whole sample obtained in both analysis is similar (see TableH)). 

5.2. B stars in Ori 0B1: Cunha et al. 

The chemical composition of the Orion star forming region de- 
rived by CL94 has been a commonly used reference for sev- 
eral applications in contemporary astrophysics for the past two 
decades. They derived LTE and non-LTE abundances for a sam- 



ple of 18 B-type stars. Both analyses were based on a photomet- 
ric determination of the effective temperature. LTE abundances 
were obtained using ATLAS6 line blanketed model atmospheres 
(iKurucd 19791) . Non-LTE abundances were derived using a sim- 
ilar strategy as in the present work, as described in Sect. 2] i.e 
Detail and Surface NLTE computations on top of a line blan- 
keted LTE atmospheric structure. Details on the non-LTE com- 
pu tations used by CL94 can be f ound i n the original papers 
bv iBecker & Butled ([T988bl [T9891 [T990allbl) and lEber & Butieil 
(Il988h . Two main differences distinguish our NLTE computa- 
tions from those used by CL94: on the one hand, older versions 
of the codes and mostly older atomic data (our model atoms are 
listed in Table [TJ, except for the Si model atom that is the same 
for both works; and on the other hand, their non-LTE lin e for- 
mation was computed on top of the slightly blanketed LTE lGoldl 
(1984) atmospheric structures, while we use the fully blanketed 
ATLAS9 model. 

In CL94, the idea of self-enrichment via supernovae ex- 
plosions was very attractive to explain discrepant abundances 
between star members of the region. This explanation to the 
chemical patterns found by CL92 and CL94 has been chal- 
lenged by new results from lSimon-Dfaz et al.l (l2006h and Paper I. 
There, highly homogeneous O and Si abundance was found for 
Ori OBI, with no clear indications of SN-II contamination. The 
high degree of homogeneity is confirmed here by means of an 
independent spectroscopic analysis. 

A comparison of C, N, and Fe abundances derived from our 
analysis, with those provided by CL94 (both NLT^ and LTE) 
is presented in Table |4] (mean abundances and 1-cr dispersions) 
and Fig. |2] (star-by-star abundances). Despite the differences in 
the methods, codes and observed spectra, the mean values of the 
present work and those from CL94 for C, N, O, and Fe, both 
in non-LTE and in LTE, are in agreement within the uncertain- 
ties. Our chemical abundances are, however, more precise thanks 
to the improvements in reducing several syste matic effects in 
the whole procedure (iNieva & Przvbillall2010allbl) . In particular. 
Fig. |2] shows significant differences in the derived abundances 
for individual stars, e.g. of up to 0.3 dex in some stars for the 
case of C and N from CL94, in contrast to our results. It is also 
remarkable the case of iron, where we find a certainly smaller 
scatter compared to their LTE determinations. Only one element 



* The NLTE abundances have been the preferred cited values in the 
literature 
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Symbols: ■ This Work (NLTE) □ Paper I (NLTE) O O CL94 (LTE, NLTE) ACHL06 (NLTE) 
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Fig. 2. Elemental abundances in our Ori OBI sample stars and a comparison with previous work, and with common abundance 
standards. Symbols for data from this work. Paper I, CL94 and CHL06 are explained in the upper part of the figure. The grey-shaded 
region represents the 1-cr dispersion of abundances from the whole sample of 13 stars. Protosolar values of AGSS09 and the Cosmic 
Abundance Standard of PNB08 are given on the right hand side of the figure for comparison. 
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has drastically changed its mean abundance: silicon. The dis- 
crepancy of Si abundance is probably related to the combined 
effect of an inaccurate determination of the effective temperature 
and the use of certain diagnostic lines, traditionally used to es- 
tablish the microturbulence and derive the Si abundance, which 
have been identified to be problematic (see Paper I for more de- 
tails). 

In a more recent work by CHL06, non-LTE Ne abundances 
for a sub-sample of stars in CL94 were derived using the previ- 
ously determined stellar parameters and s ynthetic spectra c om- 
puted with Tlusty (Hubeny 1988; Hub env & Lanzl Il995h . In 
brief, the code calculates the stellar structure and the line forma- 
tion simultaneously in NLTE. The non-LTE atmospheric struc- 
tures from Tlusty are in agreement to better than 1% with those 
from ATLAS9, in LTE, for temperatures up to (at least) ~35,000 
K (see Nieva & Przybilla 2007). Therefore the main difference 
to our work is the input atomic data for the line formation. 
Moreover, their abundances were derived from Ne i, while ours 
were determined from Ne i and Ne ii, that allowed us to analyse 
the hotter objects. A good agreement (star-by-star, and therefore 
mean abundances and scatter) is obtained between our Ne abun- 
dances and those derived in CHL06. 



5.3. B-type stars in the Solar neighbourhood 

The chemical abundances of the star sample agree very well with 
the recently proposed Present-Day Cosmic Abundance Standard 
of the solar neighbourhood, obtained from a representative group 
of early B-type stars in associations and the field by PNB08 (see 
Table |4] and Fig. |2]i. This agreement strengthens the idea that 
chemical homogeneity is naturally found in the solar neighbour- 
hood once most prominent systematic uncertainties of stellar pa- 
rameters and abundances have been minimized. These results are 
also confirmed by the study of a larger sample of B-type stars in 
the solar neighbourhood (Nieva & Przybilla, in prep.). 



5.4. The Sun 

Table m and Fig. |2] show the overall agreement between present- 
day metal abundances in the Orion star-forming region as 
derived here from early B-type stars and solaiQ values by 
lAsplund et al.l (l2009l) . despite some small remaining discrepan- 
cies. 

Comparisons of chemical composition of young early B- 
type stars and a much older star like our Sun have a histor- 
ical character, since the Sun is the preferred reference for the 
chemical abundances in most astrophysical applications. There 
are, however, several issues related to the Sun in context of 
Galactochemical evolution that are still not well understood, 
e.g. whether the Sun has actually mig rated from its nursery 
to the current location in the Galaxy (iPortegies ZwartI l2009t 
ISchonrich & Binnevll2009h . or some small-scale peculiarities in 
the solar abundance pattern in comparison with solar analogs 
dMelendez et al.1 120091: fRamfrez et al.ll2009h . The metal enrich- 
ment of the interstellar matter since the Sun's formation is also 
not precisely known, but a good estimate may be achieved via 
Galactic chemical evolution models. The agreement of chemical 
abundances derived from early B-type stars and the protosolar 
values indicate that there has been no significant local chemical 
enrichment since the formation of the Sun or that the Sun has 



been f ormed in another region - see also c.f. ICarigi & Peimberj 
(I2010h . 

We leave this discussion for the upcoming work on a larger 
sample of early B-type stars in the solar neighbourood dis- 
tributed in various OB associations and the local field (Nieva 
& Przybilla, in prep.) 

5.5. Late-type stars in OriOBI 

Late-type stars are other abundance indicators than the Sun to 
compare with our present-day abundances derived for the Orion 
star-forming region. The main differences of these stars with our 
Sun is that they have been indeed formed from the same mate- 
rial than the early B-type stars studied here and that they should 
be younger than the Sun because they belong to a young as- 
sociation. There are however a few issues to bear in mind when 
evaluating chemical abundances of cooler stars. Elemental abun- 
dances determined for late-type stars are usually not absolute 
values, like in this work, but they are derived relative to a ref- 
erence, which in most cases is the Sun. At present, atmospheric 
codes accounting for 3D effects, relevant for stars with convec- 
tive envelopes, or a non-LTE treatment of chemical elements are 
not implemented for most analyses of late-type stars other than 
the Sun. In spite of these issues, for completeness and to have a 
general picture of the Orion star-forming region, we cite a few 
recent works on the chemical composition of late-type stars in 
Ori OB I and compare them with our results. 

'Gonzalez Hernand ez et aII(l2008l) . lD'Orazi et al.l(l2009l) . and 
iBiazzo et al.. (,201ft) have recently obtained Fe abundance^ in 
late-type stars in Ori OB lb and the Orion Nebula Cluster (ONC 
- OriObld). The three studies remark the small scatter in 
the derived abundances w ithin each subgro up. For Ori OB lb, 
iGonzalez Hernandez et al.l and iBiazzo et alj derive 7 .48 + 0.09 
and 7.46 + 0.05 dex, respectivelv. while iD'Orazi et al.l find hints 
(only 1 star analysed ) for subsolar metallicity. F or stars in the 
ONC, iD'Orazi et all obtain 7. 49 + 0.07 dex , and IBiazzo et all 
7.41 +0.11 dex. As argued bv IBiazzo et al.L these discrepancies 
may be due to systematic effects in the analysis. F or example, it 
is remarkable that the re-anal ysis by Biazzo et al.l of some of the 
stars in the ONC analysed bv iP'Orazi et al.l led to differences up 
to 0. 13 dex in some cases. Our present-day abundance of iron is 
in agreement with results from the analysis of late-type stars in 
Ori OBI. 



6. Summary and conclusions 

We determined the stellar parameters and H, He, C, N, O, 
Ne, Mg, Si and Fe abundances of a sample of 13 early B- 
type stars in the OriOBI star-forming region by means of a 
detailed, self-consistent spectroscopic analysis. We applied an 
updated spectral synthesis in non-LTE based on a hybrid non- 
LTE approach using Atlas9 line blanketed LTE model atmo- 
spheres, non-LTE Detail+Surface line formation calculations, 
and line-fitting analysis techniques. We used lines from different 
ionization stages of these elements where possible, and treated 
carefully several sources of systematic errors - summarized in 
iNieva & Przvbillal (l2010albh - in order to obtain accurate and 
reliable results and minimize as much as possible the final un- 
certainties. 



' Note that protosolar abundances have been chosen because they 
represent the chemical composition at the Sun-formation time. 



The original papers provide Fe abundances on a [Fe/H] scale. We 
used the solar Fe abundances indicated in each study to compute the Fe 
abundances in the same scale used in this work. 
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The Stellar parameters derived here agree with those deter- 
mined in Paper I of this series using the non-LTE, line-blanketed 
stellar atmosphere and wind code Fastwind. The low dispersion 
of O and Si abundances in B-type stars in the Orion OBI as- 
sociation found in Paper 1 is confirmed. This is the first sys- 
tematic comparison of results derived from the codes Fastwind 
and Atlas+Detail+Surface - commonly used in the analysis 
of massive stars - from a high quality set of stellar spectra. 
Moreover, this work constitutes the first step for further com- 
parisons of elements to be incorporated in Fastwind, that is also 
able to analyse more massive stars with winds. 

The present work finds a significantly smaller abundance 
scatter (mainly in C, N, O, and Si, but also Fe) in Ori OBI early 
B-type stars than previous works and similar scatter than CHL06 
for Ne. The derived metal abundances are in excellent agreement 
with results of early B-type stars in the solar neighbourhood that 
define the Present-Day Cosmic Abundance Standard (PNB08). 
They also agree with data from an extended sample of stars in 
the solar vicinity (Nieva & Przybilla, in prep.). Despite the dif- 
ferent state of these massive young stars and of older stars like 
the Sun in terms of Galactic chemical evolution, good agreement 
with the solar values is found - which poses the question of the 
origin of the Sun. Late-type stars in the region also present sim- 
ilar iron abundances. 

This study, together with Paper I, has allowed us to establish 
a precise, homogeneous and reliable set of abundances for C, N, 
O, Ne, Mg, Si and Fe in the early B-type stars of Orion OBI. In 
a first application, the present results have facilitated constraints 
on the dust-phase composition of the Orion H ii region derived 
in Paper II. 
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Appendix A: An example of a global fit: HD 35299 

To illustrate the global quality achieved in the spectral synthe- 
sis when a final self-consistent solution is considered, we show 
in Figs. IA.3HA.7I the observed spectrum of an example star, 
HD 35299 and a synthetic spectrum computed using the stellar 
parameters and abundances indicated in Table |3] The spectral 
analysis was based on all Balmer and He lines and the metal 
lines marked at the bottom of the panels. For completeness pur- 
poses, all available lines from our present model atoms - even 
the lines excluded from the analysis for different reasons men- 
tioned in Sect. |4]- and also from other atoms not analysed here 
in NLTE are included. 

Surface gravities have been determined from the wings of the 
Balmer lines and via multiple metal ionization equilibria, that 
can be even more sensitive to this parameter. The cores of the 
hydrogen lines - which are not always well fitted - do not af- 
fect the analysis performed in this work, or at least, the effect is 
negligible. 

The line identification is based on the lines included in our 
input for Surface. Evidently, some remaining lines/elements still 
need to be included to achieve complete coverage by the global 
spectrum synthesis. Identifications of the missing lines in the 
most important blue spe ctral region may be achieved on the basis 
of the spectral atlases of iKilian et al.l (Il99lh . 
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Fig. A.3. Comparison between the observed (black) and the model spectrum (red line) for the star HD 35299, for parameters and 
elemental abundances as given in Table 3. The spectral analysis was based on all Balmer and He lines and the metal lines marked 
at the bottom of the panels. Additional lines have been included in the global spectrum for visualization purposes, but the linelist is 
not complete. See the text for details. 
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